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ABSTRACT 

The role of higher order Lyman series photons on the 21cm absorption or emission 
signature of the IGM against the background CMB during the Epoch of Reionisation 
is examined. It is shown that, taking into account the diminishing resonance line 
scattering cross section with increasing Lyman order, a non-negligible net scattering 
rate of higher order Lyman photons is expected. The resulting radiative cascades will 
substantially enhance the number density of Lya photons near a radiation source. 
It is also shown that the higher order Lyman series photons are able to collisionally 
heat the IGM by amounts of tens to hundreds of degrees kelvin. The possibility that 
the Wouthuysen-Field effect may be suppressed by the presence of dust near a galaxy 
is discussed, and it is shown that the higher order Lyman series photons would still 
induce the effect, but with a somewhat reduced 21cm radiation efficiency. It is also 
demonstrated that extended low surface brightness emission line halos will be produced 
from radiative cascades following the scattering of higher order Lyman series photons. 
These halos would provide a unique means of confirming that reionisation source 
candidates were surrounded by an IGM that was still largely neutral on large scales. 

Key words: atomic processess - cosmology: theory - line: formation - radiative 
transfer - radio lines: general - scattering 



1 INTRODUCTION 

> 

The detection of the first sources of light in the Universe 
^ ■ through their induced intergalactic H i 21cm signatu re be- 
Uh ' fore the Epoch of Reionisation (EoR) Jllogan fe Rees|[l979l : 
Pi ', IScott fc Reel Il990t iMadau et all 1 19971 ) may be realised in 
the near to not very distant future with the advent of a 
new generation of metre- wavelength scale radio telescopes, 
such as the LOw Frequency Array (LOFARj]], the Murchi- 
son Widefield Array (MWA)B the Primeval Structure Tele- 
scope/21 Centimeter Array (PaST/21CMA) 0, the Preci- 
sion Array to Probe EoR (PAPER jj, and a possible Square 
Kilometre Array (S KAfl Reviews o f this rapidly growing 
area a re provided bv lFan et al.l (|2006l ) and iFurlanetto et al.l 
(2006). Central to the detection is the decoupling of the 
spin temperature of the neutral hydrogen from that of the 
Cosmic Microwave Background (CMB). Three possible de- 
coupling mechanisms exist: coupling the energy levels to a 
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nearby bright radio source (|Bahcall fc Ekersjll969l ), estab- 
lishing thermal equilibrium with the gas through col lisions 
by other hydrogen atoms and electrons (|Fieldl TT958) . and 
coupling the energy levels to Lyman resonance line radiation 
through the Wo uthuysen-Field (W-F) effect |Wouthuvsenl 
ll952l ;fFicld 195§j), which will normally quickly establish ther- 
mal equilibrium between the spin state and the scattering 
gas. At typical intergalactic densities at redshifts z < 17, 
where the signal may be detected by the newly developed or 
planned radio facilities, coupling through the scattering of 
Lyman resonance line photons is expected to be the domi- 
nant mechanism. 

Whilst the discussion of the W-F mechanism has nor- 
mally been confined to the context of hydrogen Lya photons, 
the contribution of higher order Ly man resonance line pho- 
tons has received recent attention fearkana fc Loebl 120051 : 
lHiratall2006l : iPritchard fc Furlanettdhood h Noting that the 
higher order photons will scatter only a few times before the 
excited atom decays through an alternative channel, it was 
concluded that direct collisions by higher order Lyman se- 
ries photons would always be negligible compared with the 
Lya collision rate of photons emitted directly by the source. 
At the same time, the Lya photons produced in radiative 
cascades following the scattering of higher order Lyman se- 
ries photons would provide a substantial boost to the overall 
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scattering rate. The amount of the boost would depend on 
the maximum upper principal quantum number n (n > 3) 
0, of the photons that are able to reach a given distance 
from the source without redshifting into the the resonance 
frequency of the next lower Lyman order. 

In this paper, it is shown that the increasing mean free 
path through the IGM of higher order Lyman resonance 
line photons will result in a total scattering rate of Ly-n 
photons a few percent of that of Lya for sufficiently large 
upper principal quantum number n, much exceeding previ- 
ous estimates. The enhancement of the Lya scattering rate 
by Lya photons produced in radiative cascades is found to 
boost the Lya scattering rate by up to 30 per cent, more 
than previous estimates. It is also shown that since the 
higher order photons will not scatter sufficiently to estab- 
lish thermal equilibrium with the gas, these photons may 
provide a non-negligible contribution to the heating of the 
IGM through collisional heating. The heating may be suffi- 
cient to raise the IGM temperature in the vicinity of a source 
above the CMB temperature, resulting in an emission sig- 
nature against the CMB rather than absorption. Another 
consequence of the scattering of higher order Lyman series 
photons is the production of secondary emission lines, such 
as the Balmer and Paschen series, in radiative cascades in 
the neutral hydrogen surrounding the source. The emission 
lines would appear in the infra-red, and the emitting regions 
would subtend large angles on the sky, arising in spatially 
extended regions. Although the lines are individually weak, 
cross-correlations of the measurements of the expected lines 
may provide a detectable signal. The discovery of the emis- 
sion line halos would provide a means of confirming that 
high redshift reionisation candidate sources were in fact em- 
bedded in an IGM that was still largely neutral on large 
scales. 

Unless stated otherwise, a flat cosmology is assumed in 
this paper with a total mass density ratio to the Einstein- 
deSitter density of Q m ~ 0.3, a vacuum energy contribution 
Q v — 0.7, a baryon density Qbh 2 = 0.02, and a Hubble 
constant of H = 100/ikms" 1 Mpc" 1 with h = 0.7. 



2 THE SCATTERING OF INTERGALACTIC 
LYMAN RESONANCE LINE PHOTONS 

The optical depth through a homogeneous and isotropic ex- 
panding IGM of a Ly-n photon (with upper state principal 
quantum number n) emitted by a source at redshift zs and 
received at redshift z at frequency v > i>i u , where vi u is the 
resonance line frequency, is 



f S A > dl ? I '\ ( 1+Z ' 

Tv = o n \ dz —m(z )ipv { a„,v— — 
J z dz' \ 1 + 2 



(1) 



where ni(z') is the number density of scattering atoms 
in the lower level at epoch z' , a n — (ne 2 /m e c)fi u ~ 
0.02643/; u cm 2 Hz is the total resonance line cross section, 
fi u is the upward oscillator strength, ipv(a n , v) is the Voigt 
line profile normalized to J dvtpv(cin, v) = 1, a n is the ratio 
of the decay rate to the Doppler width Avo = Pi u b/c, where 



6 Ly/3 photons cannot decay to Lya, but will decay to Ha fol- 
lowed by two-photon emission. 



b = (2/cbT /ran} 1 is the Doppler parameter of the gas at 
temperature T and c is the speed of light, and l p is the proper 
path lengthQ In the Lorentz wing, expressed as a function 
of the normalised frequency offset x = (y — vi u )/ Avd, the 
dimensionless Voigt profile <j>v{an,x) = (Ai/n)<pv(a>n, v) is 
well-approximated by cj>v(a n ,x) ~ a„/(ivx 2 ). The differ- 
ential proper line element evolves according to dl p /dz — 
c/[H(z){l + z)] ~ (c/H )fl^ 1/2 (l + z)~ 5/2 in a flat universe 
at redshifts for which f2 m (l + z) 3 dominates the contribution 
from the vacuum energy, where H(z) is the Hubble parame- 
ter at redshift z. In the limit of scattering in the blue wing, 
the optical depth along the path of a photon emitted at 
frequency x e from a source at redshift z s and received at 
frequency x (provided it has not passed through any reso- 
nance line centre en route) is given at large separations by 



: xi 



where 



(J n a n \i u ni{z) 



(2) 



■kH(z) 

( Furlan etto fc Pritchardl l200rj : Ifflggins fc Meiksinl [2009), 
where Xi u is the wavelength of the resonance transition. 
In terms of the atomic constants, xi may be expressed 

as 



xi = 



1 / 7ve 2 \ 2 gi_ 



,2 n-i(z) 

J lu~ 



'■H{z)b' 



(3) 



where m e and e are the mass and charge of the electron. For 
large n transitions, fi u ~ (2 8 /3e 4 )n -3 ~ 1.56n -3 , so that 
xi decreases rapidly with increasing order as xi ~ n~ 6 . The 
effective mean free path of the photon as defined through 
T x = l p /lmfp, where l p (zs,z) = J zs dz'dl p /dz' is the proper 
distance travelled by the photon, is given by 
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x, (4) 



so that the resonance line photon mean free path increases 
like n 6 for large n. For sufficiently large n, the mean free path 
will match the proper distance l p for even small values of x ~ 
0(1). Thus, whilst the scattering of lower order Lyman reso- 
nance line photons occurs in the Lorentz wing, at sufficiently 
high orders the scattering will shift to the line core. The scat- 
tering of higher order photons, however, is confined increas- 
ingly to the vicinity of the source with increasing order. This 
is because, in an expanding IGM, higher order Ly-n photons 
may redshift into the next lower order. The redshift z n to 
which a Ly-n photon may travel from a source at redshift 
z s is given by 1 + z„ = (1 + 2s) [1 - (n - 1)~ 2 ]/[1 - n' 2 ]. 

The Lyman resonance order at which the transition 
to core scattering occurs may be estimated as follows. 
The Voigt line profile is well-approximated in the core by 



7 For low Lyman orders, the optical depth is complicated by 
the presence of deuterium, which will produce a Gunn-Peterson 
trough across the Ly— n profile s for x < Sl.Skms - 1 /b. For a 
D/H ratio of 2.8 ± 0.2 X 10" 5 l lO'Meara et alj|2006t) . the deu- 
terium optical depths are r a ~ 16.6 ± 0.5, rp ~ 2.6 ±0.1 and 
t 7 ~ 0.92 ± 0.03 at z = 8. The IGM becomes optically thin to 
deuterium Ly-n photons for n>3atz = 8(n>5atz = 20). 
Since the results discussed in this paper rely predominantly on 
higher orders, as the direct collisional rates of these lower order 
Ly-ri photons (above Lya) are already strongly suppressed, the 
effect of deuterium is not included. 
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Figure 1. (Upper panels) The optical depth of photons emitted 
by a source at redshift zg = 8 and received at redshift z = 7.99 
at the frequency v near the resonance line frequency vq for Lya 
through Lyn (top to bottom, alternating solid and dashed curves). 
The frequency is expressed as x = [y — vq)/Ai/]j, appropriate to 
each order. The left panel shows the optical depth for an IGM 
temperature of T = 10 K and the right panel for T = 100 K. 
(Lower panels) The values x\ of x at which t u = 1 for 0.001 < 
z$ — z < zs — z n for a source at redshift zg = 8, where z n is 
the limiting redshift at which the frequency of an emitted Ly-n 
photon redshifts into the next lower Lyman resonance frequency. 
The lines correspond, from top to bottom, to Lya (solid line) 
through Ly?; (dashed lines for higher Lyman orders). The left 
panel shows the values of x\ for an IGM temperature of T = 10 K 
and the right panel for T = 100 K. (Not included is the effect of 
deuterium; see text.) 



4>v{x) — 7T -1 / 2 exp(— x 2 ). For hydrogen gas in the temper- 
ature range 10 < T < 1000 K, the transition frequency 
x m at which the core and wing approximations match lies 
in the range 2.6 < x m < 3.3 for Lyman series photons 
with n < 8. A second order perturbation about x m — 2.8 
gives the convenient approximation x m ~ 2.8(1 + e), with 
e ~ [0.5987 - 0.1131(5.208 + loga)] 1/2 - 0.7738 (accurate 
to better than 0.1 per cent, for 0.001 < a n < 0.05). It fol- 
lows from Eq. <(3j that Lyman resonance line photons with 
n ^ 6, depending on the IGM temperature, will survive until 
redshifted into the resonance line core. In Figure [TJ the in- 
tergalactic optical depths through an IGM at temperatures 
T = 10 K and 100 K are shown for resonance line photons 
with n = 2 — 8 emitted by a source at zs = 8, along with 
the corresponding values of x\ as a function of redshift dif- 
ference from the source. It is apparent that photons with 
n > 5 or 6 will first scatter in the resonance line core. 

The source function describing the injection of Lyman 
series photons following the initial scattering of continuum 
photons by the IGM is given by 



Table 1. Probabilities p n n i for a Ly-n photon to convert into 
a Ly-n' photon per scatter. The final column gives the survival 
probability p B urv, n = Pn,n for a Ly-n photon to survive a scatter- 
ing event as a Ly-n photon. 
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nica n / dv' R n (v' ,v)n™ 
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ijHiggins fc Mei ksin 2009) ,where n™ c,n is the specific number 
density of incident continuum photons near the resonance 
line frequency v n . The incident photon number density near 
Ly-n is given by n^ nc ' n = (A-k/c)L v exp(— T v )/[(4-Kr L ) 2 hpi/] 
for a source of specific luminosity L„ and at a luminosity 
distance from the source tl ■ It is clear that the source func- 
tions of Ly-n photons depend on the order, particularly as 
a much smaller fraction of high order photons are scattered 
out of the line of sight from the source than low order. As a 
consequence, the scattering rate of higher order Lyman pho- 
tons exceeds the rate of lower order Lyman photons (above 
Lya) for a nearly flat spectrum source. 

The case of Lya photons is special, since these pho- 
tons cannot degrade. Instead, because the IGM is optically 
thick to Lya photons the photons become trapped in the 
IGM until they redshift through the Lya resonance as a 
consequence of cosmological expansion. A balance is then 
achieved between the injection rate and the escape rate of 
the photons from the IGM. The number of scatters AQcat 
before the photons escape is given by the inverse of the 
Sobolev expansion parameter 7, which is the inverse of the 
Lya optical depth in the expanding IGM, A4cat = 7 _1 = 
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Tc = X a a a ni( z)/H(z) <|Fieldlll959al ;l Gunn fc PetersorJ[l965l ; 
iMadau et aIlll997l )PI The rate of photon scatters per hydro- 
gen atom, given by integrating Eq. (5) over frequency, is then 
n^ 1 J duS^(u) = n]^ c ' a (0)ca a (a a /Tv)x^ 1 , approximating 
the Voigt profile by (a a /iv)x~ 2 . Multiplying by the number 
of scatters A/" acat gives for the scattering rate per hydrogen 
atom, P* ilcct = n^ c ' a (0)ca a since (a„/jr)(T«/j;i) = 1. Thus 
the scattering rate is simply given by the scattering rate of 
the incident photons assuming no scattering losses en rou te 
|Fieldlll959al ; I Madau et al. I I 19971 : iHiggins fc Meiksmll2009l ). 

By contrast, Lyman resonance line photons above Lya 
that first scatter in the wings will produce a negligible scat- 
tering rate because they will quickly degrade into a lower 
energy photon before they can redshift or diffuse into the 
resonance line core. A ty pical photon will not survive m ore 
than about five scatters (|Pritchard fc Furranettoll2006r i. As 
will now be shown, higher order Lyman series photons that 
first scatter in the resonance line core will produce a non- 
negligible scattering rate, contrary to previous estimates 
which assumed a constant s ource function for all t he Ly- 
man resonance line photons dBarkana fc Loebll2005l ; Ifflratal 
2006; iPritchard fc Furlanettdl2006l ). leading to a scattering 
rate estimate of order 0(5 x 10 ) times smaller than the 
Lya rate. 

The scattering rate per atom will depend on the total 
number density of resonance line photons built up in the re- 
gion in which the incident radiation field generates them. For 
a mean scattering time of 4cat,w(f) = l/nica n (j>v(a n , v) and 
a photon survival probability p SU rv,n, describing the fraction 
of Ly-n photons that produce a subsequent Ly-n photon 
upon scattering, the contribution of direct Ly-n scatterings 
from the source to the photon number density, including the 
subsequent rescatterings of the photons, is given by 



~\~ Til C&nPsurv ,ntscat ,n (^) 
roc 

x / du' R n (u', u)n*, c,n 
Jo 

+ ... 

+ (n i c<7 n p surVi „) m+1 t S cat,n(^) 

/ dv m R„(u m , u)t sca , tl „(u m ) . 



x 



dV\ R„(fi, V2)t aC a.i,,n(vi) 

du R n (u , vi)n u , 



(6) 



Here, R n (u' ,u) is the frequency redistribution function for 
the scattered photons. The survival probabilites p S urv,n fol- 



8 The equivalence .A/scat 



tering problem of Higgins 



7 also follow s from the IGM scat- 
who obtain an equi- 
librium photon number density n v = i/(njC(T a 7) in the red wing 
produced by the redshifting of photons following their injection 
by a Dirac <5- function source of strength s at a frequency fi n j, 
describing the trapping and rescattering of photons in the blue 
wing by the IGM. In equilibrium, the photons in the red wing 
must redshift away at the rate J c(i/n„/(A/scatfscat) = s, giving 
A4cat = 7 _1 - 



Table 2. The ratio Af v „ /A/",J," c (0) of Ly-n photon occupation num- 
ber to the occupation number from the incident radiation field 
assuming an optically thin IGM, divided by the IGM suppression 
factor iS„, followed by the value of the suppression factor, for a 
source at redshift zg = 8 and a T = 10 K IGM, at (proper) dis- 
tances from the source of 20 kpc (cols 2—3) and 100 kpc (cols 5—6). 
For Lya, only A/"a/A/^ nc (0) is shown. Also provided are the cor- 
responding light temperatures (T n )n due to the direct radiation 



field n 



direct, n 



(cols 4 and 7). 
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0. 


.004817 


-0.0548 








20 


4.42363 


0. 


.005624 


-0.0558 








21 


4.41521 





.006516 


-0.0569 








22 


4.40791 





.007497 


-0.0579 








23 


4.40153 


0. 


.008572 


-0.0590 








24 


4.39592 


0. 


.009745 


-0.0601 








25 


4.39095 





.011020 


-0.0612 








26 


4.38652 





.012402 


-0.0623 








27 


4.38256 


0. 


.013894 


-0.0634 








28 


4.37900 


0. 


.015501 


-0.0645 








29 


4.37579 





.017227 


-0.0657 








30 


4.37288 





.019077 


-0.0668 








31 


4.37024 


0. 


.021055 


-0.0680 









low from the spontaneous transition rates A(n,l;n' ,1') ac- 
cording to 



U(n,l;l,0), 



(7) 



where A n j = E™~ 1 1 E^.~ A(n, I; j, fc) denotes the total decay 
rate from level (n, I) to all lower levels. Values for n — 2 to 
31 are provided in Table Q] 

The direct Ly-n scattering rate is given by inte- 
grating nf,' TCCt ' n over the Ly-n cross section. Noting that 
J dvR n (y\ v) — tpv(a„, u'), the integrals in Eq. (|6]) contract, 
resulting in the direct scattering rate 



pdir 



du n v ca n ipv(an, u) 

/(I Psurv.n); 
= PA nC (0)5„/(l — p surv ,n), 



(8) 



where P; nc (0) = nt nc ' n (0)ccr n = a n L Un / ' (4nrr%hpv n ) is the 
scattering rate for an optically thin IGM, and 
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dv (fiv(a n , v) exp(-7v) 



(9) 



describes the scattering suppression factor due to the scat- 
tering out of sourc e photons by the intervening IGM out 
to th e distance tl (|Diikstra et al.l [20081 ; iHiggins fc Meiksinl 
2009). Values of S n are tabulated in Table[2]for some typical 
situations. 

Although the expression in Eq. JS} follows from a seem- 
ingly straightforward account of photon scatters, the expres- 
sion for n^ lrcct,n in Eq. ((6} reveals that an implicit non-trivial 
assumption has been made: it has been assumed that the 
photons rescattered from frequency Uj to fj+i have sufficient 
time to do so before their frequencies are Doppler shifted rel- 
ative to the expanding (or contracting gas). Since the pho- 
tons are initially scattered blueward of the line centre, red- 
shifting will generally carry the photons to frequencies with 
larger scattering optical depths. If the redshifting timescale 
(characterised by the inverse of the divergence of the veloc- 
ity field), is shorter than the scattering time at the frequency 
Vj , the redshifting timescale should be used in Eq. ([SJ . In the 
case of a contracting region, blueshifting may carry the pho- 
tons to regions of a longer mean free path, and hence longer 
scattering time. If the scattering time exceeds the contrac- 
tion timescale, then scattering may not occur at all before 
the photons are blueshifted away, and so the contraction 
timescale should be used. These factors are further compli- 
cated by the redistribution of the photon frequencies by the 
scattering, which tends to drive photons towards the line 
centre for scatters not too far out in the wing, and to con- 
centrate the photon frequencies at the incoming f requency 
for scatters well into the wings fee. iMihalasl (Il978l )). In the 
event the scattering time is shorter than the characteristic 
expansion or contraction timescale, Eq. (|6} may be used, 
from which Eq.(|8} follows. 

Incident resonance line Ly-n photons that do not rescat- 
ter as Ly-n photons will produce lower order Lyman series 
photons through the ensuing radiative cascades, adding to 
the direct scattering rate above. The cascade rates may be 
estimated similarly to the above, except that the redistri- 
bution function must now describe the conversion of a Ly-n 
photon to Ly-n' (n' < n). The Lyman photon conversion 
probabilities per scatter p n> „/ may be computed from the 
probability for an atom in state (n, I) to produce a Ly-n' 
photon through radiative cascades, 



-i j-i 



p casc (n,Z;n') = A n j ^ ^A(n,l;j, 



(10) 



The values for p casc (j,k;n') may be computed iteratively 
from low j to high, initiated by p C asc(n', 1; n') = p SU rv,n' and 
Pease (n', k; n') = for k 7^ 1. The Lyman photon conversion 
probabilities are then given by p ni „' = p ca sc(n, 1; n'). Values 
for n = 2 to 31 and n' < 5 are provided in Table [1] 

The photon density n^ lrcct '™ serves to source lower or- 
der photons through scatterings in which the original Ly-n 
photon degrades into lower energy photons. Since every or- 
der will produce further lower order Lyman photons through 
scatterings, a production cascade of Lyman photons results. 
The rates may be computed by starting with the highest 
order Lyman resonance line photons, n max , reaching a given 
distance from the central source. The production rate of the 



next lower allowed Lyman order photons is 



S. 



"maxt^n 



, (v) = nicr nm ^p 



n max ,Ti n 



(11) 



/ dv'R 
10 



n 111.1X . nmax 



/ / \ dir. 

-2{v ,v)n v , 



where R n i tn (y',v) is the redistribution function describing 
the scattering of a Ly-n' photon at frequency v' to a Ly-n 
photon at frequency v. The redistribution function is dis- 
cussed in the Appendix. The resulting density of cascade- 
produced n ma x-2(^) photons, including their rescatterings, 
is then given by 

cascade, rom^x — 2 . / \ 

n v = "iCCrn max Pn max ,n max _ 2 Iscat,n max -2(^) 

roc 

I l 1 n / I \ direct, n,-n;, 

X I dv rtn max ,n max _ 2 (^ ,V)n v , 

JO 

+ n ! CCT nmax _2Psurv,n max _ 2 

i 

/ dV\ i?n max -2(^l, V)t, 



dv R 



><:\\ x —i(y\ ) 

direct. n max 



n max ,n n 



/ / \ dir 

-i[y ,i>i)rv 



+ ... 

+ (ni C(7 n max — 2psurv,n max _ 2 J 

) 

dv m Rn max -2(Vm, // )^scat,n max _ 2 (^m) ■ 

) 

dV\ #n max -2(^l, ^2)iscat,n max -2(^l) 



Jo 



dv' R 



n max ,n n 



/ / \ din 



(12) 



As in the case of n^ lrcct,n in Eq. ©, it has been im- 
plicitly assumed that the photons are not much redshifted 
(or blueshifted) relative to the scattering medium if it is ex- 
panding (or contracting) before being scattered. If not, then 
integrating n £ ascad *> n max-s over t ne scattering cross-section 



contracts Eq. (112[) into 



pen 



direct — cascade 
,i ax -2 



Pn 



-2 ^direct 



(13) 



surv ,n n 



A similar expression results for the scattering of any directly 
generated Ly-n' photon into a Ly-n photon, 



direct — cascade 



pdr 

n' ,n 



Pn' 



I— Pa 



- Pi 



(14) 



Each cascade generated Ly-n photon in turn serves as 
a source of further lower order Lyman photons. Considera- 
tions similar to the above result in the net cascade-generated 
collision rate 



1 



(15) 



n' — n+1 



where the total scattering rate of Ly-n' photons is P n i = 

pdircct + pcascade cascade scattering rates may be 

solved for iteratively starting from higher order to low with 
P % „ = P£Z? = PC„/(1 -Psurv,n max ). For Lya pho- 
tons (n = 2), the same expression may be used except that 
1/(1 — Psurv.n) must be replaced by AT sca t = r a to account 
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for the accumulation of Lya photons as they grow in density 
once trapped in the IGM, ultimately redshifting away. It is 
the set of total rates P n that drives collisional photon heat- 
ing and the W-F effect extended to higher Lyman orders. 

It will be convenient below to express the scattering 
rates in terms of the mean photon occupation number of the 
resonance line photons averaged over the resonance line scat- 
tering profile Mv n = (c/87r)A^ J °° dvn v (j)v(cLn, v), where n v 
is the specific photon number density including the enhance- 
ments due to cascades from any higher order Lyman series 
photon scatterings and rescatterings. The scattering rate 
per atom in the lower state is then P n = {g u /gi)A u iM Vn , 
where g u and gi are the statistical weights of the upper and 
lower states, respectively, and A u i is the spontaneous decay 
rate for the transition. In terms of the occupation numbers, 
Eq. (|H| and Eq. (|15[) may be used to express the contribu- 
tions from direct scatters and subsequent cascades to the 
photon occupation numbers (for n > 2) as 



(16) 



and 



5^ A(n',l;l,0) _ 

E A( n 1-1 



Psurv.n /— ' A(n, 1;1,0) 
n — n+1 



where M„ n , = M^, cct + A^ cadc is the total occupation 
number of Ly-n' photons. For Lya photons, 



M« 



E 



A(n,l;l,0) 
A(2,l; 1,0) 



Mv n 



(18) 



The mean photon occupation number for the Ly-n transi- 
tion will also sometimes be indicated by Miu where greater 
specification of the energy levels involved is required. In par- 
ticular, for Lya photons in statistical equilibrium with the 
gas, M a i = Mti exp(— Tat/Ti), where M S i and Mu describe 
the Lya photon occupation numbers corresponding to the 
hyperfine n = 2 states i and the hyperfine ground state sin- 
glet and triplet states, respectively. For higher levels (n > 2), 
Mai = Mu may be taken. 

The spontaneous decay ra t es A(n , I; n' ,1') are computed 
following ICondon fc Shortlevl l|l97d) . Excellent ag r eemen t 
is found with the rates published in IWiese et al. fl966). 
The survival probabilities of Ly-n photons computed agree 
precis ely with those tabulated by IPritchard fc Furlanettol 
(2006), and the Lya production pro babilitie s p n ,2 agree 
precisely with those t abula ted by iHiratal ([2006) and 
IPritchard fc Furlanettol (|2006l ). who cite the fraction of de- 
graded Ly-n photons that produce Lya photons, expressed 
in terms of the values given here by p n ,2/(l — p SU rv,2). 

Representative photon occupation numbers including 
rescatterings and cascades are shown in Table [2] The val- 
ues are normalised by the occupation number of incident 
Ly-n photons assuming an optically thin IGM. A source 
with constant L\ is assumed, as this approximates a star- 
burst spectrum at waveleng ths just longward of the Lyman 
edge l|Leitherer et al.l Il999l ) . Whilst rescatterings increase 
the numbers of higher order Ly-n photons by the factor 
1/(1 — Psurv.n) above the incident number, cascades add lit- 
tle more. The exception is for Lya photons, for which cas- 
cades substantially boost the photon density because the 
Lya photons become trapped in the IGM, leaving only as 
they redshift through the resonance line frequency. Approx- 



imately one third more Lya photons are obtained compared 
with previous findings jBarkana fc Loeb 2005; Hiratal [20061 ; 
IPritchard fc Furlanett cfl2006l ; IChuzhov fc Zhend 120071 ). In- 
stead of Eq. (|18|l . these earlier estimates used the equiva- 
lent of Mr caAe = En=r PnA^/ry a ) 2 MZ(0)/(l-p snrv ,u), 
where AC" c (0) is the occupation number of the incident Ly-n 
photons in an optically thin IGM. For a source at zs = 8 in 
a T — 10 K IGM, Eq. (fT8)l gives 29 per cent, more cascade- 
produced Lya photons at a proper separation of 1 Mpc from 
the source, 32 per cent, more at 100 kpc, and 33 per cent, 
more at 20 kpc. 

It is emphasised that radiative transfer effects for Lya 
photons involving spatial diffusion, which may enhance the 
densi t y of Lya photons nea r a source (lLoeb fc Rvbickil 
1 19991 ; IChuzhov fc Zhend 120071 ; ISemelin et al.ll2007l) . or re- 
coils, which may suppress t he density near line center and 
so reduce the scattering rate (|Chen fc Miralda-Escudell2004l ; 
Iffiggins fc Meiksini r2009) have not been included. It is pre- 
sumed that the Lya photons produced in cascades will be 
subject to the same radiative transfer processes as the di- 
rectly incident photons, so that the enhancements here re- 
flect the fractional increase due to cascades alone. This is 
not obviously correct, particularly as the injected photons 
may arrive within the line core, as opposed to the wings 
as assumed in most previous studies. The actual degree of 
enhancement will depend on the relative effects of spatial 
diffusion, recoil, the local rate of expansion or contraction 
of the scattering medium, and escape from substructures 
within the IGM. 



3 COLLISIONAL HEATING BY LYMAN 
RESONANCE LINE PHOTONS 

When Lya photons are first incident on cold neutral hy- 
drogen gas, they provide a source of heat resulting from 
the mom entum transferred to the atoms by the scattering 
photons (|Madau et al.|[l997l ). This "recoil heating" persists 
for only a short period before the radiation field reaches 
statistical equilibrium with the gas through multiple scat- 
ters, and establishes thermal equilibrium with the gas (|Fieldl 
ll959d ; lMadau et al ll997l ; lMeiksinT2 006). Thereafter at most 
a residual amount of energy transfer remains in a cosmologi- 
cal setting, scaling like 7 as photons redshift through th e res- 
onance line frequency l|Chen fc Miralda-Escudel |2004| ). For 
typical IGM conditions at high redshifts, the time to achieve 
thermal equilibrium is about 1-1 yrs, correspo nding to sev- 
eral tens to hundreds of scatters (|Meiksinll2006l ); Lya photon 
collisional heating subsequently becomes an inefficient heat- 
ing mechanism. 

By contrast, because higher order Lyman resonance line 
photons do not survive long enough to establish thermal 
equilibrium with the scattering medium before degrading 
into lower energy photons, they may provide a significant 
source of heating provided their scattering rate is sufficiently 
high. The analysis for Ly-n photons that rescatter as Ly-n 
photons is identical to that for Lya photons. In addition to 
order-preserving scatterings (Ly-n to Ly-n), scatterings in 
which the incident Lyman photon degrades into lower en- 
ergy photons (eg, Ly7 into Paa, Ha and Lya), also provide 
a source of recoil heating. In this case only the photon pro- 
duced in the scattering event provides a recoil to the atom, 
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not the subsequent decays, as may be demonstrated as fol- 
lows. 

The 4-momentum of an incoming atom of rest mass 
m a and quantum state energy e, is p ai = ["/i(m a + 
ti/c 2 )vi, 7i(m a c + £i/c)]. The 4-momentum of the incom- 
ing photon which it scatters is p 7 i = (/ii//c)(n', 1). Af- 
ter the scattering event, the 4-momenta of the atom and 
photon are p a j = [jf(m a + e//c 2 )v/,7/(m c + e//c)] and 
p 1 f = (hv/c)(p., 1), for a final atomic quantum state of en- 



ergy ey. Here, 7^/ = (1 



J i,f 



'1/2 



To lowest order in 



v/c, and neglecting and e/ relative to the rest mass en- 
ergy of the atom, the resulting frequency shift between the 
outgoing and incoming photons becomes 




hvv' 



(1 



£/ - U 



(19) 



The expression reflects the change in frequency between the 
incoming and outgoing photons, including the change in the 
energy state of the atom, the Doppler shift due to the mo- 
tion of the atom, and the recoil, which depends on the ratio 
hvv' /m a c 2 and corresponds on average to an energy loss 
from the radiation field. This latter terms gives rise to the 
heating of the IGM. 

The frequency of an outgoing photon following the 
spontaneous decay of the atom after a scattering event may 
be expressed as the limit of an incoming photon with v' = 0. 
The frequency of the outgoing photon is then 

h(l — Vi • n/c) 

which expresses the energy difference of the atom and the 
Doppler shift due to its motion, without any recoil term. 
As a consequence, whilst recoil due to the scattering of a 
photon which degrades upon scattering must be accounted 
for, no recoils result from any subsequent decays of the atom 
in the non-relativistic limit. 

It is shown in the Appendix that the contribution from 
scatterings in which the incoming photon is degraded is iden- 
tical in form to that of a surviving photon: it is only the 
absorption line profile and the incident and scattered pho- 
ton energies that enter into the average amount of energy 
exchanged with the gas. (Any subsequent Ly-n photons pro- 
duced following scattering events in which the original pho- 
ton is degraded, of course, will also contribute to the recoil 
heating term if scattered, at a rate depending on the corre- 
sponding absorption line profile.) 

The resulting total heating rate due to the scattering of 
Ly-n photons is 



hv n 
m a c 2 



En — 1 
n'=l 



A(n, l;n',0) + A(n, l;n',2) 



hv„ 



A n ,i 
T 



(Tnn')n 



(21) 



for an IGM of temperature T and a harmonic mean light 
temperature of 



(Tnn')'H 



v nn > Jo" dvn v ifi V (a n ,v) 
v n / °° dv ^L-n v Lp v (a n ,v)' 



(22) 




0.01 



V 


T = 


100K 




T = 


10K 











(Mpc) 



Figure 2. Lyman photon collisional heating rate for a source 
at zg = 8 and IGM temperatures T = 100 K (solid line) and 
T = 10 K (dashed line), as a function of proper distance from the 
source. The heating rate is normalised to a Lya scattering rate 
equal to the thermalisation rate P t h ~ 6.8 X 10 — 12 [(1 + z)/9] s _1 
required to bring the spin temperature to the kinetic temperature 
of the IGM through the W-F mechanism. Also shown (dotted 
line) is the total Lya scattering rate P a /P^ c (0) normalised by 
the scattering rate of the incident Lye? photons not including the 
contribution from cascades from higher Lyman orders. 



where 



T n (v) = - — 




(23) 



fcf lMeiksinl (|2006t )1. Here v nn , = v L (l/n' 2 -1/n 2 ) (and v n = 
v„i), where vl is the frequency of the Lyman edge. A n ,i is 
the total decay rate of the p-state with principal quantum 
number n. (Note A(n, 1; n', 2) is undefined for n' < 3, and 
should be regarded as zero.) 

The frequency distribution of the photons is determined 
both by the incident radiation field and the photons pro- 
duced through rescatterings and cascades. It is convenient 
to consider the contributions from these secondary photons 
separately. 

Computing the light temperature for the secondary 
photons requires integrations over the redistribution func- 
tions. An estimate may be made most simply for Ly-n pho- 
tons rescattered into Ly-n. For scatters in the blue wing, 
the frequency redistribution may be approximated by co- 
herent scattering for which R n (v' , v) ~ ipv(a n , v)8o(v — v'), 
where 8d is the Dirac 5-function. Eq. (((J then contracts to 

n d ir oct,n ^ n mc,ft /(l _ psurv n ) = n |T'™(0)<S n /(l -Psurv.n). 

For a negligible contribution from the core, the light tem- 
perature is then given by 



r r \*]l\'l\ hVL 

2k B 



i-4 



-Xl 
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exp(-xi/x m ) 



(1 + x 1 /x m ) exp(—xi/x m ) 



-0.0338 1 



T 1/2 Xl K, 



(24) 



where x\ » x m has been assumed in the last line. The 
light temperature is independent of the gas temperature for 
wing scatters since xi ~ J 1-1 / 2 . The skewness of the photon 
number density towards the blue results in a negative light 
temperature, corresponding to a net transfer of heat to the 
gas resulting from the Doppler shifting of the photons by 
the gas. This is in addition to the recoil heating term, and 
will contribute a comparable amount. Because of the small 
value of the suppression factor S„, however, the scattering 
rate will be small, and the heating rate as well. 

A larger source of heating will arise from Ly-n 
photons that first scatter in the core x\ < x m . In 
this case, the redistribution function may be approxi- 
mated in the limit of complete redistribution, R n (p' ,v) ~ 
(pv(a n , v^tpvi&n, Eq. (JU then contracts to n dlroct,n ~ 
n 1 " c ' n +n 1 " c ' n (0)5„p S u r v, n / (1— Psurv.n) for a nearly flat source 
spectrum n^'^iO) over the line profile. The light tempera- 
ture is then 



(J-n ;h 



fc B I" 



dv exp(— Tv) 



difvifln, v) 



dv 



0.0677 1 



T 



1/2 



dx exp(- 



'(a n ,x 



Psurv,n 
1 



(25) 



dx 



K. 



Since dipy jdv < blueward of the line centre, (T™ re )H < 0, 
and the radiation field again results in a heating term in 
addition to the recoil term. The term can be large and even 
dominate the recoil heating term. Typical temperatures are 
provided in Table [2] The light temperature scales like T 1 ^ 2 . 

The contribution to the heating from photons pro- 
duced in cascades depends entirely on integrations over 
the frequency redistribution functions for the photon prod- 
ucts following the scattering of a Ly-n photon which is 
degraded into lower energy photons (see the Appendix). 
For example, in the limit of complete redistribution, 
R n / n (u' ,l>) — <pv(a n i, v')<pv(a n , v), the contribution of di- 
rect Ly-n' scatters to the number density of Ly-n photons 
produced through cascades takes the form nt 1Tact ~ cascado ' n ~ 
YZ™=l+2 n '° C ' n (°)( cr n'/ cr '0'5n'Pn',n/(l-Psurv,n). Consider- 
ations similar to the above show that the light tempera- 
ture will be very large in magnitude for a source spectrum 
n mc,n j-q^j near iy g a t across the Ly-n' line profile, so that 
these photons will heat the gas primarily through recoils. 

The net heating rate G^ at by higher order Lyman pho- 
tons is shown in Figure [5]as a function of the distance from a 
source at zs — 8 in an IGM with temperatures of T = 10 K 
and T — 100 K. The steps correspond to the increasing 
number of Lyman series transitions that contribute to the 
overall rate through radiative cascades as the source is ap- 
proached. Almost all the heating results from direct Ly-n 
photon scatterings; the contribution from cascades is about 
three orders of magnitude smaller. The heating rates are 
normalised by the local Lya photon th ermalisation scatter - 
ing rate P th ~ 6.8x 10" 12 [(l + 2)/9] s _1 l|Madau et al.lll997h . 
required to match the scattering rate of CMB photons and 
so couple the spin temperature to the kinetic temperature of 



the neutral hydrogen through the W-F mechanism. (So ge- 
ometric dilution is not included in the figure: if P a — Pth at 
a particular radius, the heating rate at smaller radii would 
be larger in proportion to 1/r 2 .) When the Lya intensity is 
sufficiently strong to initiate the W-F effect, the collisional 
heating by the scattering of higher order Ly-n photons pro- 
vides a substantial increase to the temperature of the IGM 
near a source. An increase in temperature results in a reduc- 
tion in the Ly-n IGM optical depths (see Figure [TJ , and a 
resulting increase in the heating rate. For comparison, also 
shown in Figure [2] (dotted line) is the radial trend of the 
enhanced total Lya scattering rate P a , including the effect 
of cascades from higher orders. The additional contribution 
due to cascades above the incident Lya photon rate is found 
to scale nearly like r _3//7 , although the dependence is some- 
what sensitive to the incident source spectrum. It is very 
insensitive to the IGM temperature. 

A fully accurate treatment of the heating rate by higher 
order Lyman series photons requires solving the radiative 
transfer equation for both the directly incident photons, 
with their rescatters, and the photons produced through 
subsequent radiative cascades. It is expected, however, that 
the heating rates will be reduced by at most a few per cent., 
as the photons establish only partial thermal eq uilibrium 
with the gas after only a few to several scatters l|Meiksinl 
2006), by which time they will degrade to lower energy pho- 
tons. Even if the full contribution of photons produced in 
radiative cascades is excluded, the rate due to the incident 
higher order Ly-n photons is still substantial. 



4 THE WOUTHUYSEN-FIELD EFFECT FOR 
THE LYMAN SERIES 

In this section, the Wouthuysen-Field effect is extended to 
the scattering of higher order Lyman resonance line photons. 
Computing the scattering rates at the hyperfine level intro- 
duces several complicating factors over Eq. (|6]) and Eq. (|12[) . 
First, the redistribution functions defined at the hyperfine 
level should be used. These are not straightforward, as the 
ground level is no longer sharp (in particular, the decay from 
the triplet to singlet states give rise to the 21cm line). Sec- 
ond, a photon emitted at one hyperfine resonance may red- 
shift into a redder hyperfine resonance before rescattering, as 
(Avms/v^c ~ O(0.1 kms- 1 ), requiring a travel distance of 
only 0.1 kpc at z — 8. Third, for a similar reason the absorp- 
tion profiles will greatly overlap when the thermal motions of 
the atoms are included. Fourth, the survival probabilities of 
individual hyperfine resonance line photons must be used. 
Fifth, scatters in the wings complicate the averaging over 
the Voigt profiles since the natural line broadenings of indi- 
vidual hyperfine transitions differ from the total line widths 
for a Ly-n transition. This latter effect moreover results in 
a light temperature in the context of the W-F effect that 
differs from the light temperature of Eq. (|22|l in the context 
of photon collisional heating, unli ke the case fo r scattering 
rates dominated by core scatters (lMeiksinll2006l ). 

To make an estimate of the magnitude of the role higher 
order Ly-n photons play in the W-F mechanism, these ef- 
fects will be neglected. In particular, it will be assumed that 
the photons enter the core so that the Voigt profile shapes 
for the hyperfine transitions are the same as for their par- 
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ent Ly-n transitions. Since the scattering rates of photons 
that scatter far in the wings are negligible, this should be a 
reasonably good approximation to the overall effect. Accord- 
ingly, the photon occupation numbers for a given hyperfine 
transition become the same as for the parent Ly-n transi- 
tion. The hyperfine scattering rate per atom in the lower 
state is then P„ = (g u / 'gi)A u iNu n , where g u and gi are the 
statistical weights of the upper and lower hyperfine struc- 
ture states, respectively, and A u i is the spontaneous decay 
rate for the hyperfine transition. 

Because the level populations above the ground state 
are assumed to be excited exclusively by the scattering of 
Lyman resonance line photons, in equilibrium the ratio of 
the occupation number of the excited states to the ground 
state levels will all be of the order of the photon occupation 
numbers of the incident Lyman photons or less, depending 
on the branching ratios from the upper state. The 2 2 S state 
is an exception. Since radiative decays from the 2 2 S state 
to the ground state are forbidden by the selection rules for 
dipole transitions, this state will fill up until other processes 
become effective. The dominant mechanism for vacating the 
state is through two-photon decays. For this state, the oc- 
cupation number will be larger than the other excited states 
by a factor on the order of the ratio of the dipole decay 
rate to the two-photon decay rate, still resulting in occu- 
pation numbers much smaller than the ground state levels 
because of the smallness of the photon occupation number. 
Accordingly, resonance line scatterings from states above the 
ground state may be neglected, as these will result in cor- 
rections to the state populations quadratic in the photon 
occupation numbers of the incident radiation field. Denot- 
ing the singlet and triplet state occupations by n s and nt, 
respectively, and all the hyperfine levels above the ground 
state by 7ij with i = 1 to i = N indicating the levels that 
may be reached by scattering Lyman resonance line photons 
up to order n, the radiative cascade equations become 



dn s 
~dT 

dn t 
~~dt 

and 



;=i 

N 



i = l 
N 



= n ^ ~ n tYl irMiAit + Ps%s - P£n t (26) 



= giAisMsiUs + -giAitNurit + ^ rijAji — 

i=i+i 

ru i^Ais + An + Y A ilJ > ( 27 ) 

where M s i and Afu denote the photon occupation number 
between state i and the singlet and triplet ground state, re- 
spectively, Aij denotes the electric dipole Einstein A sponta- 
neous decay rate from hyperfine level i to hyperfine level j, 
and gi — 2Ft + 1 denotes the degeneracy of hyperfine level i, 
where Fi — Ji + 1 for total angular momentum Ji — Li + S, 
nuclear spin / = ±1/2, angular momentum Li = l iy where 
U is the orbital angular momentum of state i, and electron 
spin S — ±1/2. The transition rates are computed similarly 
to the transitions above, but applying the Russell- Saunders 
multiplet formalism twice to allow for J and I coupling to F. 
The resulting rates were checked by verifying that applying 



the sum rules recovers the fine-structure and full sponta- 
neous decay rates. 

It is convenient to label the states by increasing n, then 
by increasing L followed by increasing J and F. Since only 
spontaneous decays from levels n to levels n' with n' < n are 
considered, this labelling ensures the transitions from states 
with higher labels to lower will always correspond to tran- 
sitions from higher energy states to lower. (It is noted that 
the labelling scheme does not correspond always to increas- 
ing energy levels. For example, the n 2S+1 Lj = 2 2 S\/ 2 en " 
ergy level is hi g her th an both the hyperfine levels of 2 2 P 1 / 2 , 
iBreit fc Tellerl (1940)). Since fine and hyperfine transitions 
between states with the same principal quantum number are 
neglected, there is no need to reorder the labelling by energy 
within these quantum structure levels. If transitions up to 
n — n max are allowed, then the total number of hyperfine 
states above the n = 1 states is N = 2(n^ ax — 1). Labelling 
the 2 2 S°/ 2 state by i = 1 and the 2 2 S\, 2 state by i = 2, the 
rates A\ s and A 2 t are taken to correspond to 2-photon emis - 
sion at the rate A 2l — 8.23 s -1 l|Spitzer fc Greenstein|[l95l] ) . 
The corresponding two-photon absorption terms from the 
ground state singlet and triplet states to the pair of 2 2 S\/ 2 
hyperfine states are therefore highly negligible and excluded. 
Eqs. ((26} also include radiative excitation and de-excitation 
between the singlet and triplet hyperfine states, given by 
the respective rates P^ and P^, to allow for an incident 
continuum radiation field due to the CMB or nearby radio 
sources. Similar collisional terms may be added if required. 

The excited levels n > 1 will rapidly achieve statistical 
equilibrium on a timescale of order A~f ~ 0(1O -8 s) (except 
for the 2 2 5*1/2 states, for which the timescale is on the order 
of 0.1 s). The states may thus be assumed to have reached 
a steady state, so that Eq. (|27[l reduces to a matrix equa- 
tion. It is convenient to renormalise the level occupations by 
n» = ni/[n s A/}, nc (0)], where ATj^ c (0) is the incident photon 
occupation number just longward of the Lyman edge, and 
to define the ratio r — n t /n s between the triplet and singlet 
ground state levels. Denoting the rescaled occupation levels 
by the vector n, the matrix equation for the excited levels 
becomes 



Mn = y, 



(28) 



where M is an upper triangular matrix with diagonal el- 
ements Ma = A is + A it + 53i=i-^«> Mij = — Aji and y 
is an absorption vector with elements y\ = y 2 = (indi- 
cating no 2-photon absorption to states i = 1 and i — 2), 
and m = giAi s Vsi + (gi/2,)A it rvti for i > 2. Here the no- 
tation v sl = A/" si /AC c (0) and v ti = A/"ti/A^ nc (0) has been 

introduced. 

Following iFieldl (Il958l ). the value for r may be deter- 
mined by casting Eq. (|26|l in the form 

" Xy n t - Pj7n, + Pun t - P^n 3 , 



dt 
dnt 
~~dT 



= R 



= P. 



Ptl n t + P^n s 



Punt , 



(29) 



where P^ y and P^/ are the effective excitation and de- 
excitation rates of the triplet state by the scattering of Ly- 
man resonance line photons. Defining T st as the effective 
temperature corresponding to the 21-cm transition, given 
by hcj(k,B A2icm) = T st — 0.068 K, and the 21-cm spo nta- 
neous transition rate A ts = 2.85 x 10" 15 s" 1 ( Wild 1952| ), in 
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equilibrium, r = 3exp(— T a t/Ts), where Ts is the spin tem- 
perature, is given by r = (P s 



Ly + Pj?)/(P t L / + PS), where 



Pjf = 3A ts (rfl/r s t) and P£ = A ts (l + Tr/Th). Here 
is the brightness temperature of the incident continuum ra- 
diation field at the 21cm transition frequency. In statisti- 
cal equilibrium without an incident continuum (Tr = 0), 
r = 3exp(-T S i/r L ) = P s L //P t L /, where T L is defined as the 
light temperature, and is identical (to order k^Tst/hVa) to 
the harmonic mean temperature Eq. (I22p of the radiation 
field when the sca ttering rate is dominated by core scat- 
ters {Meiksin 2006|) . It relaxes to the kinetic temperature of 
the gas after a few dozen to a few hundred scatters a cross 
the resonance line centre (|Fieldlll959bl ; iMeiksinl 12006") . Us- 
ing the definition of Tl above, the spin temperature may be 
expressed generally as 



10 




%t Tst (1 + T R /T at ) (Tl + VL y 

T L ° S Tr exp(T s t /Tl ) /Tl + J/Ly 



(30) 



15 20 

redshif t 



where yi^ y = T st P ts y /{Ti,Ata)- (The expression must be 
modified if electron or atomic collisions are important, here 
assumed negligible.) 

When the light temperature of the Lya photons has 
relaxed to the kinetic temperature of the gas, the problem 
is entirely determined once the triplet de-excitation rate P^/ 
is specified. Fo r pure Lya s cattering, P t ^ y = (4/27)P a , where 
P a = 3Af a A a (|Fieldlll958l ). The higher order generalisation 
is given by 

n uF \ a / 

x A(n,l,J,F;l,0,±,l)p s (n,l,J,F), (31) 

where Aftj n> i,j,F) is the occupation number of photons with 
frequencies at the hyperfine structure resonance between 
the ground state triplet state and the state (n, I, J,F), and 
Ps(n, I, J, F) is the probability that an electron in the hyper- 
fine state (n, I, J, F) produces a decay to the ground singlet 
state, 

p 3 (n,l,J,F) = A~l JF A(n,l,J,F;n',l',j',F') 

n' <n,l' ,J' ,F' 

x p 3 (n',l',J',F'). (32) 

Here, A nUF = £„/<„,!/, j, >F , A{n, I, J, F; n',l', J',F') (= 
A n ,l by the sum rules), is the total decay rate of state 
(n, I, J, F). The probabilities p„ may be solved for iteratively 
from Eq. (32}, initiated by p s (2, 0, ±,0) = 1, p s (2,0, |, 1) = 
0, p.(2, 1, |, 0) = 0, p,(2, 1, |, 1) = 1/3, p.(2, 1, |, 1) = 2/3, 
and p s (2,l,|,2) = 0. 

The solution to the equations shows the singlet-triplet 
excitation and triplet-singlet de-excitation rates are boosted 
only by the increase in the Lya collision rate P a due to the 
Lya photons produced in cascades following the scattering 
of higher order Ly-n photons (Figure [2}. The ratio P t ^ y = 
(4/27)P Q is found to be accurate to a small fraction of one 
per cent, even close to a source (down to a proper separation 
smaller than 20 kpc for a source at zs = 8). 



Figure 3. Lyman photon collisional heating rate at a distance 
of 100 kpc from a lOAfoyr -1 starburst, for source redshifts in 
the range 6 < zg < 26. Shown for IGM temperatures of T = 
100 K (solid line) and T = 10 K (dashed line). Also shown is the 
corresponding total Lya collision rate P a /P^ nc (0), normalised by 
the collision rate of the incident photons. 

5 ASTROPHYSICAL CONSEQUENCES 

The redshift z r of the reionisation of the IGM is constrained 
primarily by measurements of intergalactic Lya absorption 
in the spectra of high redshift QSOs and by polarisation 
measurements of the CMB. The spectra of high redshift 
QSOs constrain z r > 5.7, a bove which the Lya optical be- 
comes immeasurably large (|Fan et al J 12006). The five-year 
Wilkinson Microwave Anisotropy Probe polarisation data 
yield z r = 11.0 ±1.4 assuming sudden reionisation, with 
2cr an d 3a lower limits of z r > 8.2 and z r > 6.7, respec- 
tively (|Dunklev et alj|2009h . Radio telescope efforts to de- 
tect the EoR through its 21cm signature are focussed on 
appropriate redshift ranges accordingly: 6.1 < z < 10.8 
(LOFAR), z < 17 (MWA), 6.1 < 2 < 27 (PaST/21CMA) 
and 6.1 < z < 13 (PAPER). Since the impact of higher 
order Lyman photons on the IGM is sensitive to redshift, 
the astrophysical consequences are illustrated for a range of 
redshifts. 

5.1 Photon collisional heating 

The impact Lyman photon collisional heating may have 
on the temperature of the IGM is illustrated in Figure [3] 
using a starburst galaxy. A galaxy continuously forming 
low metallicity (Z = 0.05^©) stars at the rate lOMQyr -1 
with masses between 1 < M < WOMq and a Salpeter 
IMF will have a steady luminosity at A = 915 A of L v = 
3.8 x 10 28 ergs -1 Hz" 1 after about 10 r yrs (Leitherer et all 
1999). This corresponds to a photon occupation number at 
the Lyman edge of AO — 4.8 x 10 _23 rj^p C , where the (proper) 
distance is normalised to 1 Mpc. 

The heating rate is found to decline with increasing 
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source redshift as the IGM becomes increasingly optically 
thick to Lyman photon scattering, reducing the effective Ly- 
man photon collision rates. A second reduction factor with 
redshift is the decrease in the number of Lyman orders that 
may contribute at a fixed proper distance before a photon 
redshifts into the resonance frequency of the next lower or- 
der. Each step in Figure [3] corresponds to a decrease by one 
order. For comparison, also shown in Figure [3] (dotted line) 
is the redshift trend of the enhanced total Lya scattering 
rate P a /P™ c (fi), including the contribution from cascades 
following the scattering of higher order Lyman photons. The 
rate is normalised by the incident rate excluding the cascade 
contribution. 

Because the optical depth of the IGM diminishes with 
increasing temperature, the photon collisional heating rate 
strengthens with increasing temperature, varying as G^ at oc 
T 1 / 2 . This gives rise to a weak thermal instability that grows 
quadratically with time. Photon collisional heating may thus 
produce a substantial increase in the temperature of the 
IGM surrounding a source and transform an absorption sig- 
nature against the backgr ound CMB into an emission sig- 
nature (|Madau et al.lll997f ). 

5.2 W-F effect in the absence of Lya photons 

It was found that the radiative cascades following the scat- 
tering of Ly-n photons will enhance the number density of 
Lya photons by as much as an order of magnitude near 
a source. Near a galactic source, however, dust may be 
present. It is unclear that Lya photons will be able to sur- 
vive sufficiently long to build up the photon density re- 
quired to initiate the W-F effect before being absorbed 
by dust grains. Observational searches for dust surround- 
ing galaxies suggest dust may be present at least as far 
out as 100 kp c arou nd low redshift galaxies |Zaritskvll 19941 ; 
iMenard et al.l |2009| ) . If comparable amounts of dust were 
introduced into the surroundings of galaxies at the time of 
the EoR following bursts of star formation, then the long 
total path Lya photons must travel to migrate from the 
blue frequencies at which they become trapped in the gas 
at large distances into the Doppler core, where they be- 
come effective scatterers, may excee d the mean free path for 
dust scattering. The measurement of lMenard et al. I d2009h of 
A v ^ 2.4 x 10 _3 (r p /100/i _1 kpc) -0 ' 84 at a projected separa- 
tion of r v suggests a mean free path for colliding with dust 
grains of about 4 x 10 5 (r/100ft _1 kpc) 1,74 kpc a distance r 
from the centre of the galaxy, assuming Ax — 1.086NdQe<Jd 
for a column density of Nd dust grains with cross sections 
ad and an ext inction efficiency factor of Q e ~ 2 in the visual 
l|Spitzerlll978T ). For Lya photons with a mean free path on 
the order of the distance to which they are able to survive 
before becoming trapped within the blue wing of an atom, 
the distance travelled by the photons (of order N sca , t — 1000 
in the wings times the Lya mean free path), before diffus- 
ing into the resonance line core approaches the mean free 
path for scattering off a dust grain, which may cause some 
suppression of the W-F effect. The density of dust around 
starbursts, however, is likely to be considerably larger than 
in a normal galaxy, by as much as tw o orders of magnitude , 
and may extend to quite large radii l|Heckman et al.| [2000). 
In this case the diffuse Lya photon radiation field required to 
drive the W-F effect may become substantially suppressed. 
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Figure 4. (Upper panel) The de-excitation rate of the ground 
triplet state, normalised by the rate corresponding to the thermal- 
isation scattering rate P t h, induced by higher order Lyman series 
photons in the absence of Lya photons. The rate is shown at a 
distance of 100 kpc from a 100 Mq yr _1 starburst, for source red- 
shifts in the range 6 < 25 < 26. An IGM temperature T = 10 K 
is assumed. (Lower panel) The corresponding 21cm radiation ef- 
ficiency 77. 



Only the photons that have scattered too few times to have 
migrated far into the core would contribute. 

The loss of the Lya photons, however, does not mean 
that the W-F effect will be extinguished. Although far fewer 
in number, the higher order Lyman series photons will still 
induce a W-F effect. Because the higher order photons will 
not achieve thermal equilibrium with the gas, the spin tem- 
perature will no longer reflect the kinetic temperature of the 
gas but rather the net exchange rate between the ground 
state singlet and triplet levels resulting from collisions by 
the higher order Lyman photons. 

The evolution of the triplet-singlet de-excitation rate 
at a (proper) distance of 100 kpc from a lOOM0yr _1 star- 
burst (as above), is shown in Figure U (upper panel) for 
an IGM temperature of T = 10 K. (The results are found 
to be very insensitive to temperature.) The de-excitation 
rate is only a fraction of the rate for a total Lyman pho- 
ton scattering rate matching the thermalisation rate Pth, 
but is non-negligible. The light temperature Tl is found 
to be high, varying only moderately with redshift, from 
T L = 170 K at z = 6 to T L = 200 K at z = 26. The spin 
temperature established, however, must include the scatter- 
ing of CMB photons. Shown in the lower panel is the re- 
sulting 21cm radiation efficiency factor 37 = 1 — Tcmb/Ts, 
which controls the strength of the 21cm signal through 
ST ~ ?7T S (1 - e~ T )/(l + z), where ST is the change in the 
antenna temperature compared with th e CMB temperatur e 
in a medium with 21cm optical depth r dMadau et al.|[l997l ). 
Typical efficiency factors of only one per cent, or smaller re- 
sult, so that a detection would be difficult with the upcom- 
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Figure 5. (Upper panel) The de-excitation rate of the ground 
triplet state, normalised by the rate corresponding to the ther- 
malisation scattering rate P t i 15 induced by higher order Lyman 
series photons in the absence of Lya photons. The rate is shown 
for a 100 Mq yr — 1 starburst at 25 = 8. An IGM temperature 
T = 10 K is assumed. (Lower panel) The corresponding 21cm 
radiation efficiency r). 



ing telescopes. The efficiency, however, rises at smaller radii, 
as shown in Figure [S] for the galaxy at zs — 8. It is noted 
the timescale for the Lyman resonance photons to affect the 
spin temperature will be on the order of (P t ^ y ) _1 — 1 Myr, 
so that bursts of shorter duration may not have much ef- 
fect. A brighter source, however, would increase the signal 
proportionally. Also, scatters by the residual Lya photons 
still outside the resonance line core, not included in the dis- 
cussion here, will contribute as well, further enhancing the 
signal. 

Although the emphasis in this paper is on the effect 
of higher order Lyman photons around a single source, a 
diffuse radiation field will result in regions sufficiently near 
sources that higher order photons from multiple sources are 
able to penetrate them before scattering through a lower 
order Lyman resonance. In this case the collective incident 
radiation field from the sources is given by 

Jv = \hj z dz ' ii? (1+ % {ev,{z,) exp [" r -( z ' *')]>' ( 33 ) 

l|Meiksinll2007l ). for a source emissivity e„/, where v' = ^(1 + 
z')/(l + z), and r„(z, z') is the total optical depth due to all 
the Lyman resonances (and any other contributing terms) 
encountered by the photons from gas between the redshifts z 
and 2'. The indicated spatial average is over the sources and 
the IGM optical depth jointly, as these are likely correlated 
since sources form in overdense regions. 



Table 3. Scaled emission line photon number emissivities e nn i = 
n n A nn i (in units of s — 1 ), for a source at redshift zg = 8 and a 
T = 10 K IGM, at (proper) distances from the source of 20 kpc 
(cols 2-3), 100 kpc (cols 4-5) and 1 Mpc (cols 6-7). The results are 
shown in pairs of Balmer and Paschen lines with upper principal 
quantum number n. 
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58433.7 
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10273.5 


6 


41193.9 
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32874.8 
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7984.15 
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6931.77 
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6782.75 
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18464.0 


6665.82 
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6571.86 






17 
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18065.2 


6494.71 
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18717.1 


6806.52 


17912.9 


6430.56 






19 
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6720.97 
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6649.25 
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6588.95 
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6537.39 
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18008.3 


6493.12 
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17919.5 


6454.62 
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17841.7 


6421.08 
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17772.5 


6391.43 










27 


17710.7 


6365.05 










28 


17655.4 


6341.53 










29 


17604.7 


6320.09 










30 


17558.5 


6300.60 










31 


17515.8 


6282.66 











5.3 Cascade radiation emission lines 

A consequence of the scattering of Lyman series photons is 
the generation of a host of intermediate emission lines pro- 
duced in the cascades following the scattering. The emissiv- 
ities for these transitions may be computed from the state 
occupancies given by the solution to Eq. (|28p . 

The scaled Balmer and Paschen photon number emis- 
sivities, e„2 = n n A n 2 and e„3 = h n A n 3, summed over the 
hyperfine structure lines for a given upper principal quan- 
tum number n, are shown in Table [3] for a source at zs — 8 
at (proper) distances of 20 kpc, 100 kpc and 1 Mpc from 
the source, for an IGM at temperature T = 10 K. (The 
emissivities are very insensitive to temperature. The values 
change by no more than a few per cent, for T = 100 K.) The 
bolometric emissivities in physical units are then given by 
eij — (hh / '4)huijA/Y lc {0)uj , where the factor of 4 takes into 
account that only 1/4 of the ground state hydrogen atoms 
are in the hyperfine singlet state (see §12}. The emissivities 
€ij for a source at zs = 20 are about a factor 5 smaller. 
The emissivities are found to diminish slowly with distance 
from the source. The lower order Balmer sequence emissiv- 
ities scale to better than 8 per cent, accuracy as r- _1//3 for 
Ha, r _1//4 for H/3 and r _1//s for H7. These scalings apply to 
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Figure 6. Cascade radiation emission line profiles around a 
lOOM0yr _1 starburst galaxy at zg = 8 in an IGM with 
temperature T = 10 K. (Upper panel) Shown are bolomet- 
ric Ha, H/3 and H7 intensities, corresponding to the respective 
observed wavelengths of 5.9/nri, AAfim and 3.9/xm. The inte- 
grated fluxes within one arcminute of the source are /h q = 
1.6 X 10- 18 ergcm- 2 s- 1 , / H/3 = 8.2 X 1CT 



/ H7 = 5.8x10 ia erg cm 
sity at 9 = 10 arcsec. 



(Lower panel) Ha specific inten- 



the Paschen series Paa, Pa/3 and Pa7 as well. The halo of 
any given emission line vanishes beyone the horizon of the 
lowest order Ly-n required to generate it. 

The bolometric intensity a projected distance b from a 
source at redshift zs for a bolometric emissivity varying as 
eij =e*,(r/r»)- Q - 2 is 



2tt lJ 



1 r(i±g) . 

27rV2 Q r(a/2) ^ 



<*+ 2 (l + 2 S )- 4 / d/(& 2 + / 2 )-( Q + 2 >/ 2 (34) 



f) < 1 + **>~ 



(35) 



where the redshift factor has been introduced to convert the 
bolometric intensity to the value measured at z = 0. The 
corresponding specific intensity is 



1 



47T H(z) Vij (1 + 2 S ) 3 



+ 



H{z)r* 



-(a + 2)/2 



where Aw is the velocity difference from line centre. A homo- 
geneous and isotropic medium around the source has been 
assumed. More realistic line profiles would need to include 
the structure of the underlying density and velocity fields. 

The profiles may be illustrated by the cascade line ra- 
diation produced by a IOOM0 yr" 1 starburst at zs = 8 in a 
T = 10 K IGM. The bolometric intensities for Ha, H/3 and 
H7 are then i Ha oi 1.1 x 10~ 20 9~ 4/s erg cm -2 s" 1 arcsec -2 , 
in/3 — 4.3 x 10~ 21 #~ 5//4 erg cm -2 s _1 arcsec -2 and ih 7 — 
2.6 x 10~ 21 #~ 6/ ' 5 erg cm -2 s" 1 arcsec -2 , at an observed an- 
gular separation from the source 9 measured in arcsecs. 



The profiles are shown in the upper panel of Figure [6] 
In the lower panel is shown the specific intensity at 9 = 
10 arcsec, i v ~ 2.5 x 10~ 32 erg cm" 2 s" 1 Hz" 1 arcsec" 2 [1 + 
(Aw/47 km s" 1 ) 2 ] " 7/6 . 

Although the intensities are small, the integrated fluxes 
from the emission line halos around the source are appre- 
ciable. At zs = 8, Ly/3 photons are able to travel as far 
56 Mpc from the source before redshifting into the Lya res- 
onance. This corresponds to an angular distance of 3.4°. 
The total integrated flux is ~ 5.5 x 10" 17 erg cm" 2 s" 1 . 
Whilst isolating the flux from such an extended halo is likely 
unfeasible, even within the inner arcminute the total flux 
is /h<*(1') — 1.6 x 10" 18 erg cm" 2 s" 1 . Similarly, the fluxes 
for H/3 and H7 are /h/j(1') ^ 8.2 x 10~ 19 erg cm" 2 s" 1 and 
/h-y(I') — 5.8 x 10~ 19 erg cm -2 s _1 . An observing strategy 
involving narrow band near infra-red imaging over a wide 
area would be required to detect the halos. The existence 
of multiple emission lines has the advantage that it would 
be possible to correlate the signals in frequency space, with 
continuum emission from contaminating sources subtracted. 
For example, defining Sfnp = /h/3 + /cont - (/cont), and sim- 
ilarly for Hy, the correlation of the H/3 and H7 signals is 
(SfnpSfBi) = (/H/3/H7) + [(/cont) - (/cont) 2 ], so that the 
signals only need to be detectable above the noise of the 
continuum across the frequency band of the detector. The 
emission lines will also give rise to a diffuse sky brightness. 
At a few such bright sources per square degree, a minimum 
near infra-red bolometric sky brightness of ~ 0.1 mjy sr" 1 
may be expected; the value would be much larger if a much 
higher density of sources was required to reionise the IGM. 
Although the detection of individual hyperfine emission lines 
is unlikely in the near future, it is noted that the state oc- 
cupancies at the hyperfine structure level are found to scale 
nearly in proportion to their statistical weights 2F + 1. 

It was assumed in §|4]that the optical depths of all levels 
above the ground n — 1 levels were negligible so that absorp- 
tions from them could be ignored. This assumption is now 
justified. As discussed in §|4] the ratio of the occupancies to 
the ground state occupancy will be on the order of the pho- 
ton occupation numbers, and typically much smaller. Only 
the 2 2 5i/2 states have large values, of 10 3 — 10 5 , because of 
the smallness of the two-photon decay rate. 

The optical depth for transitions between hyperfine 
states i and j is (for absorption from i to j), 



Xij aijrii(z)H(z) 
1549 - j 



(37) 



9i \ ^aA a 



iA/i nc (0)(l + z) 3/2 



(36) 



where XfjAji has been normalised by the value for Lya. 
Even for the 2 2 Si/2 hyperfine states (ni and ri2 in the no- 
tation of this paper), the optical depths will be tiny. Bal- 
ancing n\A2~, ~ N™ c (0)A a n s gives an upper limit hi ~ 
0(A a /A2j — 8 x 10 7 ); in practice the value is more than 
two orders of magnitude smaller because of the suppression 
factors S„ and reduced cross-sections for n > 2. Since gen- 
erally hi < hi, using this generous upper limit for hi shows 
that Tij < 1 provided 



\ 6 A 

X~A a 



8.48 x 10" 12 [AC c (0)] '(! + . 



-3/2 



(38) 



an equality that is readily met for any relevant transition 
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and realistic radiation field. As an illustration, for a source at 
zs = 8 in a T — 10 K IGM, at a distance of 100 kpc from the 
source the largest n ^ n hyperfine transition optical depth 
is found for absorption from (n,l,J,F) = (2, 0, i, 1) — + 
(3,1,|,2) with t ~ 1.3 x 10 10 A^ I1C (0) ~ 5 x 10" 11 for 
A/"a nc (0) ~ 3.5 x 10~ 21 corresponding to the thermalisation 
Lya scattering rate Pth at z = 8. 



6 CONCLUSIONS 

Higher order Lyman series photons from sources at high red- 
shift may have a substantial influence on the 21cm signature 
from the IGM in the vicinity of the source before the IGM is 
reionised. Taking into account the diminishing cross section 
for scattering with increasing order, it is found that Ly-n 
resonance line photons with n > 5 or 6 will first scatter in 
the resonance line core, where they scatter efficiently. The 
consequences for the subsequent rate of collisional heating, 
the intensity of the Wouthuysen-Field effect, and the forma- 
tion of extended diffuse halos of emission lines produced in 
radiative cascades are summarised. 

Whilst Lya photons readily establish thermal equilib- 
rium with the IGM, higher order resonance line photons 
degrade into lower order photons before they are able to 
do so. Because of the much smaller fraction of higher or- 
der Ly-n photons scattered out of a line of sight compared 
with Lya, they are able to produce a substantial amount 
of collisional heating. The heating arises primarily from the 
blue distortion of the photon density distribution, as only 
photons blueward of the resonance line frequency are able 
to survive to large distances without being scattered out of 
the line of sight. Typical light temperatures for Ly-n pho- 
tons, where n is the principal quantum number of the upper 
state, of (T„) H ~ -0.015T 1/2 K are obtained, where T is the 
temperature of the IGM. Resulting collisional heating rates 
of several tens to hundreds of Kelvin per Gyr are obtained. 
The heating rate scales as T 1//2 , inducing a weak thermal 
instability in the IGM, with the IGM temperature increas- 
ing quadratically with time. One consequence would be that 
warm regions are heated more rapidly than cooler regions, 
which would contribute to the patchiness of a 21cm absorp- 
tion signature against the CMB where the spin temperature 
is still less than the temperature of the CMB. 

Radiative cascades following the scattering of higher or- 
der Lyman photons will add substantially to the number 
density of Lya photons near a source, enhancing their num- 
ber by an order of magnitude. The number of Lya photons 
produced in cascades is found to be as much as 30 per cent, 
higher than previous estimates. 

Measurements of the dust content in nearby starburst 
galaxies suggest the W-F effect induced by Lya photons 
may be partially suppressed due to the absorption of Lya 
photons by dust grains, and possibly even completely elim- 
inated, in the vicinity of starburst galaxies at high redshift 
during the EoR if comparable amounts of dust were present. 
It is demonstrated that higher order Lyman series photons 
will induce a non-negligible W-F effect even in the complete 
absence of Lya photons. The 21cm radiation efficiency at 
distances exceeding 100 kpc from the galaxy would then be 
one per cent, or smaller except for extremely strong bursts, 
although the signal is stronger closer to the galaxy. The sig- 



nal would still likely be too weak to be discovered around an 
individual galaxy by first generation 21cm detection experi- 
ments, although it may be detectable by a Square Kilometre 
Array. Future radio detections and upper limits may provide 
a means of exploiting the suppression to constrain the dust 
production rates and transport efficiency in high redshift 
galaxies. 

Cascade radiation following the scattering of higher or- 
der Lyman photons will produce emission line halos extend- 
ing over several arcminutes around a source, limited only 
by the distance to which Lyman photons may travel before 
redshifting into the resonance frequency of the next lower 
order. It is shown that the integrated lower order Balmer 
fluxes within one arcminute surrounding starburst galaxies 
would be substantial. The detection of the halos would pro- 
vide a unique means of confirming that candidate reionisa- 
tion sources are in fact surrounded by an IGM that is still 
largely neutral. 

The effects discussed here are in the context of a ho- 
mogeneous and expanding universe. Accounting for cosmo- 
logical structures will m odify the results in several ways 
l|Higgins fc Meiksinl 120091 ) . The local density and velocity 
fields will affect the mean free paths of the Lyman photons, 
as will heat input from other possible sources such as Active 
Galactic Nuclei or Quasi-Stellar Objects. Of particular in- 
terest are H n regions surrounding photoionisation sources. 
The horizons of high order Lyman photons could then be 
substantially extended as the ionised gas would become op- 
tically thin at their resonance frequencies. The collisional 
heating, W-F effect, and the production of emission line 
halos would then extend into the neutral gas just beyond 
the ionisation front. Even with these complications, it is ex- 
pected that the salient features of the effects discussed here 
will still be present. 



APPENDIX A: REDISTRIBUTION FUNCTION 
FOR DEGRADED LYMAN RESONANCE LINE 
PHOTONS 

The redistribution function for an incoming Ly-n' and an 
outgoing photon other than Ly-n' differs from the normal 
redistribution function for scattering photons in which the 
outgoing photon is the same as the ingoing photon (Ly-n' — * 
Ly-n'). The principal difference arises from the re-emission 
function following absorption. In the case of a preserved Ly- 
n' photon, the outgoing frequency is the same as the ingo- 
ing in the restframe of the atom. In the case of a degraded 
photon, the outgoing freque ncy distribution is given by the 
appropriate Lorentz profile (|WeisskopJll933h . 

As an example, a scattering event is considered in which 
a Ly-n' photon is degraded into a final Ly-n photon, with 
two intermediate transitions (the minimum required to pro- 
duce a final Lyman photon) . The final photon need not be a 
Ly-n photon; the example only serves to illustrate the com- 
putations involved. 

Accordingly, four levels are considered in the atom, la- 
belled 0-3, corresponding respectively to the ground state 
(n = 1), the lowest energy excited state, a higher inter- 
mediate sta te, and the n' en ergy level. The formalism and 
notation of IWeisskopJ (|l933h are followed. The absorption 
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profile for an incoming photon of frequency £' in the frame 
of the atom is 



/(A = 77* 



(Al) 



where the line-centre absorption frequency is 1/3 and the to- 
tal decay width of the upper level is S3. The decay from 
state 3 to state 2 involves recoil. Because the upper state 
was excited from a definite energy level (the ground state), 
the upper energy level is definite. As a consequence, only the 
uncertainty in the energy level of state 2 will determine the 
width of the emission profile. Since the atom will decay to 
the centre of energy level 2 on average, the frequency of the 
emitted photons will peak at £ = [h£' — (£2 — e o)]/h — f — 1/2, 
where vi = (e2 —ta)/h. The probability density for emitting 
a photon of frequency £ in the restframe of the atom is then 
given by 



P3,2 (£',£) = - 



12 



(A2) 



where 72 is the decay width of state 2. The redistribu- 
tion function in the restframe of the atom is given by 
/(£')P3.2(£'i £)• Note that if the incoming photon is absorbed 
blueward (redward) of the line-centre, the emitted photon 
frequency £ will also be shifted blueward (redward) of the 
line centre frequency V23 = (£3 — (-2)/ h. 

For an atom moving at velocity v in the laboratory 
frame, the redistribution function in terms of the incoming 
and outgoing photon frequencies v' and v in the laboratory 
frame becomes, allowing for the Doppler shifting of the fre- 
quencies to first order in v/c, 



R v {v' , ft'; v, ft) = 



[„'_„ 3 (l + vn7c)p+^ 



(A3) 



where the recoil term Ai> re coil = (/i"3"23/w, a c 2 )(l — ft ■ ft') 
has been included. 

It is convenient to define the dimensionless frequencies 



"23 



V - "3 



with Avd = fa-, (A4) 
c 



where the frequencies are normalised by the same Doppler 
width. The dimensionless decay widths 03 = 73/AfD and 
0.2 = 72/(z / 23&/ c ) are a l so introduced, along with the recoil 
parameter e — I1V23/ m a bc and the dimensionles velocity u = 
v/6. It is also convenie nt to introduce the new coordinate 
system (|Hummerlll962r ) 

fti = 7+(n' + ft), ft2=7-(fi' — ft), 63 = ft x ft', (A5) 

where 7± = [2(1 ± /i)] -1 ^ 2 and fi = ft ■ ft'. A dimension- 
less velocity-dependent redistribution function is defined by 
R u (x', ft'; q, ft) = (Avd) 2 Rv{v' , ft'; v, ft), where q = x — x' . 
It is given by 



Ru(x', ft'; q, ft) = 



V 7T J (x 1 



{x> - 7+V/2 - l^u 2 /2) 2 + o 2 



V23C12IV3 



:7- w 2 



1-1 I, "23 



"23 
"3 



+ £(!-/*) 



+ 



"23 
"3 



-0 2 



(A6) 



It is convenient to perform the velocity average over a 
Maxwellian of the redistribution function in Fourier space. 
The Fourier transform of Eq. (|A6[) is given by 



Ru(k, ft'; A, ft) 



j / ikx 

ax e 



dq e lXq R u (x' ,n;q, ft) 



,1.-1 / "2 . 
exp^ -i7 + «i ( k — — A 



+ -«7_ u 2 



K — I 1 + — — I A 

"3 



"23 



— iAe(l — u) — asl/tl — d2 1 Al 

"3 



(A7) 



Averaging over a Maxwellian velocity distribution gives 



R(k, ft'; A, ft) 



-3/2 



d 3 ue u R u (k, ft'; A, ft) 



exp 



II 1 2 "23 I , I 
— 0,3 K ft 0,2 A 

4 1/3 



"23 



ie(l-/x)A+- K A-- l + A 



1 "23 

2 V3 



(k — X)\\i 



(A8) 



Moments of the dimensionless frequency shift are given 



by 

a n (x,fi) = 



q n R{x',n;q,h) 



(A9) 



The Fourier transforms of the moments may be expressed 
as 



S»(K,/i) = 



d{i\) n 



R(k, ft'; A, ft) 



(A10) 



noting that the factor q n corresponds to the n th derivative 
with respect to i\ in Fourier space. The transformed mo- 
ments are readily computed from Eq. (|A8|) . Angle averaging 
the results and Fourier transforming back to frequency space 
give for the first three moments, 

Oo(a;) = (/>v(a3,x), 

, . 1 d(j} V (a 3 , x) 
ai{x) = 2-^x £ ' 



1 



a 2 (x 



1 "23 



I + \ ^^(03, x) 



r(a 3 ,x) 



3 vi 



dx 2 



dx 

(All) 



where only the first order terms in e are retained. The first 
order moment ai (x) is of the sam e form as for the scattering 
of Ly-n to Ly-n (Meiksin 200(J), and results in the heating 
rate given by Eq. (|21[) . The Fourier transform of the redis- 
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tribution function averaged over angle is 
sinh ^jj 



R3,23(K, A) 



"3 V "23 



2 "3 V "23 



x exp 



i.i 1 2 1/23 i\i 

-83 /t -K 02 A 

4 V3 



■ , 1 , 1 

ie\ H — kA 

2 4 



1 + 4 1 ^ 



(A12) 



where the subscript on -R3,23(k, A) indicates the scattering 
of a photon of frequency vz to V2z ■ 

Similar redistribution functions may be defined for the 
decay product photons of frequencies v\i and V\, however 
the analysis become considerably more complicated as tran- 
sitions betwe en non-sha r p leve ls are involved. Following the 
discussion of IWeisskop J l| 19331 ). the frequency distributions 
may be multiply peaked and correlated with the frequencies 
of previously emitted photons. For example, for photons ab- 
sorbed and emitted in the Lorentz wings, the distribution 
of photons produced in the transition from level 2 to level 1 
could have two peaks. In one, the emitted frequency would 
be correlated with the frequency vp of the photon emitted 
in the transition 3^2, peaking at v' — vp — v\, followed by 
a transition 1 — > with the photon frequencies peaking at 
v\ — [e\—eo)/h. A second peak would occur at the frequency 
^12 = (fi2 — ei)/h, with the transition 1 — > producing pho- 
tons peaking in frequency at the difference between V2 and 
the frequency of the photon emitted in the 2 — > 1 transition. 
As the most probable frequency for the final emitted Lyman 
photons (in this example) is independent of the incoming fre- 
quency v' , a good approximation to the redistribution func- 
tion Ly-n' to Ly-n is R„i , n (v' , v) = (pv{a n i , v')(pv(a n , v). 
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